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Abstract 

This paper is devoted to description of the relationship among oriented associativity equations, 
symmetry consistent conjugate curvilinear coordinate nets, and the widest associated class of semi- 
Hamiltonian hydrodynamic-type systems. 
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1 Introduction 

The nonlinear partial differential system 1 

Q2 c i Q2 c m Q2 c i d 2 c m 



daJda m da k da n da k da m daWa n 



(1) 



describing a displacement vector appears in [4]. Following [14] we call this system the oriented associa- 
tivity equations. This system admits the scalar linear spectral problem (cf. [4]) 

d 2 h d 2 c m dh 



da % dai da l da^ da 

or, alternatively, the vector linear spectral problem (cf. e.g. [15] and references therein) 



da k ^da k da m ^ 

Integrable system (1) was extensively investigated in a number of papers (see, for instance, [9]) dedicated 
to the so-called coisotropic deformations. Some other aspects were considered in [8], [11], [12], [13], [15]. 

This paper is devoted to integrability of system (1). We consider a connection of (1) with a widest 
class of semi-Hamiltonian hydrodynamic-type systems 



Be 1 

(4) 



we present a geometrical interpretation for (1) and linear spectral problems (2) and (3), and describe some 
transformations preserving (1), (2) and (3). 

The celebrated WDVV equation (see, for instance, [4], [5], [6], [10]) 

d 3 F d 3 F d 3 F d 3 F 

tfrm _ ^mn j k (5) 



da' l daWa m da n da k da s da % da k da m da n daWa 

can be obtained from (1) by the potential reduction (where r] ks is a constant nondegenerate symmetric 
metric) 

In this paper, we follow [4], and step by step unravel the relationship among (1) and the widest class of 
the so-called conjugate curvilinear coordinate nets determined by (see, for instance, [1]) 

diPjk = PjiPik, i^JT^k, (7) 



5(3 lk = 0, i^k, (8) 

where d{ = d/dr l ,5 = S<9 m is the so-called shift symmetry operator, and the rotation coefficients (3 ik de- 
pend on N Riemann invariants r k . The first subset of the above equations is the famous Darboux system, 
while the second condition means that the rotation coefficients depend only on differences of the Riemann 
invariants. For the sake of simplicity we shall call these conjugate curvilinear coordinate nets symmetry 
consistent. 

^^Here and below the sum over any of repeated indices in the opposite locations (that is, one subscript and one superscript) 
is understood; otherwise the sum is indicated explicitly. 
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The paper is organized as follows. In Section 2 we introduce the metric and the basic set of solutions of 
linear systems (determining the symmetry consistent conjugate curvilinear coordinate nets) whose compat- 
ibility conditions imply the oriented associativity equations. In Section 3 we construct the transformation 
from the symmetry consistent conjugate curvilinear coordinate nets to the oriented associativity equations. 
In Section 4 we construct the relationship between linear spectral problems for the oriented associativity 
equations and the symmetry consistent conjugate curvilinear coordinate nets. In Section 5 we equip the 
oriented associativity equations by the unity condition. We concentrate here on the three-component case. 
In Section 6 we prove that these oriented associativity equations are a Hamiltonian system. In Section 7 we 
prove that the oriented associativity equations can be interpreted as a system of equations describing N- 
component position vector of a hypersurface in centroaffine geometry. In Section 8 we construct the inverse 
transformation from the oriented associativity equations to the symmetry consistent conjugate curvilinear 
coordinate nets. In Section 9 we construct an infinite set of particular solutions for semi-Hamiltonian 
hydrodynamic type systems whose rotation coefficients depend on differences of the Riemann invariants 
only. Finally, in Section 10 we just emphasize the relations among these very important systems. 



2 Linear Spectral Problems. Basic Sets of Solutions 

Consider two linear systems 

diH k = /3 ik Hi, dii) k = £ w Vi» i ^ k > ( 9 ) 

whose rotation coefficients /3 ik depend on differences of the Riemann invariants r k only (see (8)). This 
means that a particular set of solutions H k , ^ satisfies two extra equations 

5Hi = XH U dtpi = X^. (10) 

The first set of compatibility conditions dj(diH k ) = di(djH k ),dj(diip k ) = di(djif) k ) leads to a full set of 
equations (7) describing conjugate curvilinear coordinate nets, while the second compatibility conditions 
d 5 {5H k ) = 5{d 3 H k ) } d^ k ) = 6(dM yield (8). 

Remark: The symmetry consistent conjugate curvilinear coordinate nets are well known in classical 
differential geometry (see e.g. [1] and [17]). Moreover, system (7), (8) was derived more recently in the 
context of algebro-geometric solutions for multidimensional integrable systems (see [3]). It is interesting 
to note that this system also arises in quantum statistical physics, see Slavnov [16], and in theory of the 
discrete analogue of conjugate curvilinear coordinate nets known as D-invariant lattices (see [2]). 

In this case N infinite series of solutions of systems (9) can be recursively found by quadratures (see 
[17]) 

5Ht l+1 ' s) = H^ s \ 54 n+1 ' s) = 4 n > s \ , = 1,2,..., AT, n = 0,l,..., (11) 

where Sifjf'^ = and 8H^ ' S ' = 0. 

Choose N particular solutions H( k )i = H as the basic set of solutions. Then introduce a non- 
degenerate (and non-constant in the generic case) symmetric metric 

N 

9ik = y^#(i)m-^(fc)m- (12) 
m=l 

Lemma: N particular solutions ip^ = ipf'^ can be chosen in the form 

4 S) = g sn H [n)u (13) 

where g sn is an inverse metric for g~i k . 
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Proof: Substituting (13) into the second equation of (9) yields 

d^H^k) = P kl g sn H {n)ll i ± k. 

Upon removing the parentheses on the l.h.s. and multiplying both sides by g~j s the above relations boil 
down to 

H {n)k g ns dig sj = {/3 ik - l3 ki )H ij)h i ^ k. 
Taking into account (recall that SH^i = 0) the equations 

diH {s)i = -}^ mi H( s ) m , diCj jk = yXPim ~ Pmi)( H (j)i H (k)m + H(k)iH(j) m ), (14) 

one obtains an identity. The theorem is proved. 

Then, obviously, three additional identities (where 8\ and 5i k are the Kronecker symbols) 

N N N 

3*=I>SM!?, *i = X)^S%)m, 5 lk = Y,4 m) H(m)k (15) 
m=l m=l m=l 

hold. In the Egorov case (3 ik = (3 ki and the metric g lk = const (see (14)). 

3 Reconstruction of Oriented Associativity Equations 

In contrast with the previous section, we consider N Riemann invariants r k as functions of N independent 
variables t n , i.e. we introduce N commuting hydrodynamic-type systems 

r\ k = %&ri, (16) 
Hi 

where Hi is an arbitrary solution of the first linear system in (9). These hydrodynamic-type systems are 
semi-Hamiltonian (i.e. possess infinite set of conservation laws parameterized by N arbitrary functions of 
a single variable, see [17]). In such a case, they can be written in the conservative form 

d t kh = d x g k , (17) 

where dih = ipiHi, dig k = ^i-H(fc)i and ip i is an arbitrary solution (parameterized by N arbitrary functions 
of a single variable, see [17]) of the second system in (9). Introduce N conservation law densities a k such 

that di<x k = tp[ k ^Hi and N conservation law densities c k such that diC k = Hi, where Sipf* = t/j[ k ^ (i.e. 

ijf ] =4 1,k \ see (11) and (13)). 

Theorem: iV commuting hydrodynamic-type systems (16) can be written in the conservative form (4), 
whose compatibility conditions are oriented associativity equations (1). 

Proof: The relation ij)f^ = c k ^p\ s ^ follows from 

N N 

;.(*) 



dc k = Y$mH m dr m = c k da s = H m dr m , 

m=l m=l 

where c\ = dc l /da k (see (4)). Taking into account the second identity from (15), one can obtain 

4 = X^S^Wm- (1* 



m=l 



Then (see (9)) 



U) m 



If we take into account that the expression in the first brackets vanishes because SHq^ = 0, and the 



expression in the second brackets is nothing but ip\ (since 6$F=,I>\ k) ), then one can conclude that 



dic k = ^H {3)1 . (19) 
On the other hand, if the hydro dynamic- type systems (16) possess N conservation laws (4), then 



u m u i t k — <J m c k i x . 



Substituting the r.h.s. of (16) for r l k leads to (recall that dia k = ipf^Hi) 



94 = ^l^H, = 4 k) H U} 



i ■ 



which coincide with (19). The theorem is proved. 

Remark: The second derivatives of the functions c* with respect to the field variables a\ a k can be 
easily derived from (19). Indeed, the relations 

d lC ) = c k Js d ia s = r'jj/r". = 4 k) H {j)i 

lead (upon multiplying the third and fourth blocks of the above expression by the ratio H^/Hi, and 
summing according to the second formula from (15)) to 



J m n {j) m n (k)m 

m=l 

which, obviously, satisfy (1) by virtue of (15). Its symmetric form 

N 

c ijk = 9is c jk = ^ ^ 



q, = E ^Y"" . (20) 

11 1 



in 



m=l Hm 

is well known in the theory of WDVV associativity equations and Frobenius manifolds, where such expres- 
sions for c^k can be found using the theory of meromorphic differentials on algebraic Riemann surfaces 
(see, for instance, [1] and [10]). 
Remark: Since 

m=l m=l 

one can easily derive (recall (15) and di<i k = i/jf^Hj) 

(21) 



'da k ~ ~H~' 
5 



Thus, the characteristic velocities i>o%(a) = H^yjHi (see (16)) of N commuting hydrodynamic-type 
systems (4) are nothing but dr l /da k . Thus, the compatibility conditions dv^/da^ = dv l ^/da k should be 
satisfied. Indeed, multiplying both sides by af and summing over i, one obtains an identity, which follows 
from 

r)nP N rinP 

E m uu m _ V^.,m uu m 

m=l m=l 

where da l m /da k = a l ms v^ k y Thus, we conclude that 

dr l 

^ = ^)( a )- ( 22 ) 

Any semi-Hamiltonian hydrodynamic-type system (see [17]) 

a\ = vial, (23) 

is associated with the non-degenerate metric tensor g lk . The necessary and sufficient conditions for exis- 
tence of this tensor are given by (here is the covariant derivative) the Tsarev lemma (see [17]) 

g ik vi = g jk vi V,r k V,rf. (24) 

However, for hydrodynamic-type systems (4) we should not solve this system, because we already know 
that the metric tensor in the field variables a k is given by (12). Indeed, we have (see (12), (15)) 

N N N N N N / N \ 

i=l k=l i=l k=l s=l n=l \m=l / 

N N N / N \ / N \ N 

s=l n=lm=l \j=l / \fc=l / m—1 

So, we arrive at the conclusion that the diagonal metric gkk = H k in the Riemann invariants, in perfect 
agreement with Tsarev's definition of semi-Hamiltonian metric (see [17]). 

Anyway, taking into account the expression for the covariant derivative Viv\ = div\ — T^i^ + Y\ m v 
conditions (24) for hydrodynamic-type systems (4) reduce to a more compact form 



in 
k ■ 



rfkj _ -jk A -pfc m _ pfc m /«r\ 
y u ks ~ y ^ksi 1 im u js ~ 1 jm u is- ) 

By virtue of (15) and (20), the first group of equations yields an identity. To consider the second group of 
equations, at first we should compute 

r i - l 7l i™( d 9mk dg mj dg jk 



2 a V daj dak da 
Taking into account (14), (15) and (21) yields 



m=\s=\ ^ ' s=l s p=l 

Then, upon taking into account (15) and (20) once again, the second group of equations in (25) becomes 
an identity. 
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4 Linear Spectral Problems for Oriented Associativity Equa- 
tions 

As it was already mentioned in Introduction, the oriented associativity equations (1) can be obtained as 
compatibility conditions of the scalar linear spectral problem (2) 

h ik = Xc s ik h s , (26) 

where we denote c l jk = d 2 c l /da j da k , hi = dh/da\ h ik = d 2 h/da l da k . 

Lemma: The function /i(a) is a generating function of conservation law densities (see (17)). 
Proof: If h(a) is a conservation law density, then (see (21) and cf. (18)) 



r dr m N H (k)m N 



daJ ~ ^ ' ~H tn 

m=l m=l " L m=l 

Also, taking into account (26) and (19), one can compute 

dhk = hk m da m = Xc s km h s da m = Xh s dc s k 

N N N N N / N \ 

= XY^Y,^nH {s)n )(^H ik)m dn = I>«»^? ^nH ( k )m dr m , 

s=ln=lm=l n=lm=l \s=l / 

i.e. we have 

N 

dh k = \Y^i) m H {k )mdr m . 

So, we arrive at the relation 

\m=l J 

which immediately reduces to the second equation in (10). The lemma is proved. 

Now, consider the commuting hydrodynamic-type system (r is a group parameter) 

< = & a k T = bl (27) 

where Hi is an arbitrary solution of first linear system (9) and dib k = ip\ Hi. As it was mentioned in 
Introduction, oriented associativity equations (1) can be obtained as compatibility conditions of the vector 
linear spectral problem (3) 

b\ = K s b s , (28) 

where we denote b\ = db t /da k . 
Lemma: The functions 

N 



b k (a) = X-^^Hn 



m=l 



are conservation law fluxes for generating functions of commuting flows (27). 
Proof: Taking into account (28), (19) and (15), we find that 



N / N \ / N 



db k = \c k ns b s da n = Xb s dc k = Y^^™ Hm ^^n )H (s)ndr n 



s=l \m=l / \n=l 
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N N / N \ N 

= EE E^^ #^-^ n = E^- drr 

m=ln=l \s=l / m=l 



Thus, indeed, <9jfc fc = ^ On the other hand, we have 



(k) 
i 5 



i.e. 



Xdtf = H^Pim^ + (sH t - ^PrniHm ) + V^JjS^m + ^ f <ty?° " E^- ) ■ 

Since the first sum equals to the fourth sum, and the second sum equals to the third sum, we arrive at the 
conclusion that Xdib k = ip^SHi (recall that Si/)} = 0), which agrees with the previously computed dib k = 
ijj^Hi, if and only if Hi satisfies the first equation in the linear spectral problem (10). The lemma is proved. 
Remark: One commuting flow (27) can be found without integration, i.e. (y is a group parameter) 

a; = {a s c\ - c%. (29) 

Indeed (see (20)), we have 



di(a c s - c ) - a c sn dia - a c sn ip i hi-a^^ u Vi H * 

1 J- m. 



E Wm' ll {s)m I ST^ An) rj \fv s i (k) tt 
?v Zjl>\ H{n)m \H i = a >H {s)i . 



N N jMjt tt 
n=lm=l 



Since di(a s c k — c k ) = tp i Hi, where Hi is some solution of the first system in (9), we conclude that 
Hi = a s H( s )i. Substituting this expression into the first system in (9) leads to an identity, and the 
substitution into the first equation in (10) yields 



N N N / N 



6H = H {s)l 5a s = J> (s)i J>^tf m = E 5>S?ff(.)i H m = H % . 



s=l m=l m=l \s=l 



Thus, the family of hydrodynamic-type systems (4) has a simple commuting flow (29), which in the 
diagonal form (see the first equation in (27)) has characteristic velocities H/ H such that SH = Hi. 



5 Reduction to Canonical Form 

Hydrodynamic-type systems (4) have N additional natural conservation laws 

dQ 



l* = 9 X -^- (30) 



da k 



Indeed, suppose that the hydrodynamic-type systems (4) possess N additional conservation laws c\ k 
d x Q\. Then dQ\ = c l m dc™, i.e. Q\ s = c l m c™ s . Thus, we conclude that Q l ks = Q\ k and N functions Q % 
determine the fluxes of these conservation laws. The compatibility conditions dQ l ks /da? = dQ l j s /da k hold 
by virtue of (1). 
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Rewrite N commuting hydrodynamic-type systems (4) in the differential form, i.e. (see (18)), 

N N N 



k=l m=lk=l 



dz 1 = a l dx + ^V fc (a)<it fc = a l dx + ^^^ m H{k)m,dt k . 



Thus, we see that functions 4. (a) do not depend on the choice of Lame coefficients H^, while (recall again) 
did k = ifj\ k ^Hi. Below we shall consider the reduced version, i.e. 



AT N 



= Y,Y,^H {k)m dt k . 



m=l k=l 



In such a case, we pick the first "time" variable t 1 and choose the new field variables a 1 = c\(sl), and then 
instead of (4) we shall consider just N — 1 commuting flows 

* = < 31 > 

where dc l /dd k = dc l /da k , k = 2, 3, N. Thus, all functions c l can be found in quadratures, i.e. 

N 



771=1 

Indeed, rewrite the conservation laws (30) in the differential form 

N N 

dy< = c%si)dx + J2Qk(*)dt k = c\a)dx + Q\{a)dt l + ^g^.(a)dt fc . 

k=l k=2 

Choose c % = Q\(a), then dc l = Q\ s da s = c l m c™ s da s = c % m dc!^ = c % m dd m . So, dc l /dd k = dd/da h . 
Obviously, the hydrodynamic-type systems (31) can be written in the diagonal form 

^ = f^r f Y (32) 
#(i)< 

Indeed, the first of these commuting flows (16) can be written in the form 

i Hi i 
r t i. 



X TT 

Then all other commuting flows (16) reduce to (32) upon substituting for r l x from the above formula. 

Below we shall omit tildes over field variables a k and c n , because we are going to consider just the 
oriented associativity equations supplemented with the so-called "unity" condition (cf. [4]) 

4 = 4, (33) 

which is equivalent to considering a canonical set of iV — 1 commuting flows (31). We shall call such 
oriented associativity equations normalized. 

Thus, the conservative representations (31) reduce to the form 

a] k = d t ml, a\ k = d t i (a 1 ^. + u\), i,k = 2, 3, N, 
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where the new unknown functions u n (a , a , ...,a ) appear from integration of (33), i.e. 

c 1 = ^(a 1 ) 2 + d' = a 1 a k + u k , k = 2,3,...,N. (34) 

In this case, the shift symmetry operator S can be easily expressed via the field variables a k instead of the 
Riemann invariants r n , i.e. 

N d N ( N da k \ d N ( N \ d d 

5 = E^ = E E^T 9^ = ElE^W™ = ^T- (35) 

m=l fc=l \m=l / k=l \m=l / 

Thus, two equations of scalar and vector linear spectral problems (26) and (28), namely 

h x = Xh, b\ = Xb k (36) 

coincide with the eigenvalue problem for the above shift symmetry operator, while the remaining equations 
in (26) and (28) become, respectively, 

N 

h lk = \J2<khs, i, k = 2, 3, . . . , N, (37) 

s=l 

b{ = X ( (51 - 5[8\)b l + f2 u Lb m ) , * = l,2,...,iV, fc = 2,3,...,JV. 

\ m=2 / 

Example: Consider the simplest nontrivial case N = 3. Then two commuting flows are given by (here 
x = t 1 ,t = t 2 ,y = t 3 ) 

at = 9 x u b , b t = d x (a + v b ), c t = d x w b ; 

(38) 

a y = d x u c , by = d x v c , c y = d x (a + w c ), 

where a = a 1 ,b = a 2 ,c = a 3 ,u = m 1 (6, c),f = u 2 (b,c),w = u 3 (b,c) and the subscripts indicate the 
corresponding partial derivatives. The compatibility conditions for (38), which read (a y ) t = {(h) y, {b y )t = 
(b t ) y , (c y )t = (c t ) y , lead to three algebraic equations relating the second-order derivatives 

Ubb = v bc Wbb - v h bW bc + w 2 bc - w bb w CCJ u bc = v cc w h b - v bc Wbc, u cc = v\ c - v bb v cc + v cc w bc - v bc w cc , 

which are integrable by the inverse spectral transform. The scalar Lax pair for the latter (see (37)) is 
given by 

/ o i o \ / h a \ / h a \ / o o l 

= A ^66 h , hb = A U bc V bc Wbc 

b \ u bc v bc w bc ) \ h c ) \ h c J c \ u cc v cc w cc 
while the adjoint Lax pair has transposed matrices, i.e. 

u bc u cc 

A | V bc Vcc 

1 W bc W cc 

Remark: Both commuting hydro dynamic-type systems (37) also have the same Lax matrices, i.e. 

u bc u cc 

v bc v cc 

1 Wbc W cc 











Ubb 


Ubc 








Vbb 


Vbc 








Wbb 


Wbc 







{» 


b 2 


V 1 


b 2 


b 3 ) 




V b 3 










Ubb 


Ubc 


i 


i-i 


i 


Vbb 


Vbc 








w bb 


Wbc 









)l 


i 
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6 Hamiltonian Structure of Oriented Associativity Equations 

The system of quadratic equations 

U b b = v bc w bh - v hb w bc + w 2 bc - w bb w cc , u bc = v cc w bb - %%, u cc = v\ c - v bb v cc + v cc w bc - v bc w cc , 

is nothing but the oriented associativity equations in three-dimensional case with the unity condition (33) 
(see also (34)). Introduce a new set of field variables q 1 = u bb , q 2 = u bc , q 3 = v bb , q 4 = v bc , q 5 = w bb , q 6 = w bc . 
Then this quadratic system becomes a six-component hydro dynamic-type system 

12 2 o <f<f + <?V ~ Q 2 Q 3 

q c = ?6» i c = d b - 5 , 



3 4 4 ~ (f + Q 4( l 6 

Qc = Qb, Qc=d b - 5 , (39) 



5 6 6 o (g 6 ) 2 - <?¥ + ? 4 9 5 - g 1 
Qc = Qb> Qc= d b ~ 5 • 

As it was already mentioned in Section 5, in the present paper we restrict ourselves to the generic case 
when all characteristic velocities v l ^(a) are pairwise distinct. We seek the Riemann invariants r(a,b,c) 
for both commuting hydrodynamic-type systems (38) written in the diagonal form 

r t = V(2)r x , r y = v {3) r x . (40) 

Taking into account that r a = 1 (thanks to the shift symmetry operator 5 = d a , see (35)), one can obtain 

r b = v {2 ), r c = V( 3 ), (41) 

where the characteristic velocities vm and vt$\ of (38) are related polynomially, 

(v { 2)) 2 - q\ 2 ) ~ Q 1 
%) = ~ 5 • 

while V(2) satisfies the characteristic equation (45) for the first hydrodynamic-type system from (38), i.e. 

(^ (2) ) 3 - (q 3 + q 6 )(v {2) ) 2 + (gV - q 4 q 5 - q 1 ^ + q\ & - q 2 q 5 = 0. 

Since the three characteristic velocities are distinct, the commuting hydrodynamic-type systems (38) can 
be written in diagonal form (40), where the Riemann invariants can be found by quadratures (see (41)), 
i.e. (cf. (46)) 

r k = a + j \v\ 2 )db + v 1 fe ) dc), k = 1,2, 3. 

Remark: Hydrodynamic-type system (39) possesses at least two additional local conservation laws. 
Indeed, the existence of three Riemann invariants (see (41)) implies three additional conservation laws 

<9 C W( 2) = d b ^^ -g-^ , k = l,2, 3. 

Hence, the roots v%\ are conservation law densities. However, five conservation law densities q 3 , q 6 , v^,v 2 2 ^, 
related by the linear equation 

q 3 + q 6 = vf 2) + vf 2) + vf 2) (42) 
11 



by virtue of the Viete theorem. This means that just two of the above three conservation laws are new. 
Main Result of this Section: Upon expressing q l ,q 2 and q 6 using the Viete theorem via q 3 ,q A ,q 5 

7 (2)' V (2)i V (2)> 



and vL},v%s,v%-., the hydrodynamic-type system (39) can be written in the local Hamiltonian form (cf. [6]) 



'ds l 



where the flat coordinates are 



S ~ V {2)i S - 


- «(2), 


s 3 = 


U(2) 




s 4 


= d\ 


the metric is 
















( - l 
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1 








\ 
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-1 


1 











_ 1 
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1 


-1 

























2 
















2 
























1 / 


the momentum density P 


is 












i 1- 


i k 












q = s = 




+ s 2 


+ 


s 3 


) 2 -( 



2g 3 -(^2) + 4) + 4)^ 



( 1 1 \ 
10 10 
1 1 
1 
1 
\000002/ 



and the Hamiltonian density H is 



sV + sV + sV - - a 

4 



6 ^ 2 - sV 



9 2s 1 s 2 s 3 + (s 1 + s 2 + s 3 - s 6 V 
q= 2? ' 



while the other conservation law densities q 3 , g 4 , q 5 , g 6 (recall (42)) are just linear combinations of the flat 
coordinates s k , i.e. 



q 3 = -(s 1 + s 2 + s 3 + s £ 
2 



Cf = s \ g 5 = s 5 ? g 6 = I (s l + s 2 + s 3_ s e 

2 



7 Centroaffine Geometry 

Thanks to the presence of the shift symmetry operator (i.e. hi = Xh, see the first equation in (36)) linear 
spectral problem (37) for normalized oriented associativity equations (see (33)) becomes 

N 

h ik = >^u s ik h s + X 2 u\ k h, i, k = 2, N. 

s=2 

This linear system arises in centroaffine geometry. Following [5], consider a linear overdetermined 
system (which is required to be compatible of rank N) 



(43) 



for the ^-component position vector r = r(a 2 , . . . , a ) of a hypersurface M N ~ 1 in centroaffine geometry, 
where A is a spectral parameter, gij(a 2 , . . . , a ) is a pseudo-Riemannian metric, and T^(a 2 , . . . , a ) are 
components of a torsionless affine connection (which in general is not the Levi-Civita connection of flat 
metric g^). The conformal class of gij is nothing but the second fundamental form of M N_1 . 
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The compatibility conditions for this system have the form (see [5]) 



F s ij9sk-T s jk gsi-T s ik g sjj d k gij-dig jk - djg ikj - diT s jk - djT s ik , 

(44) 

WTk + (JuV = r° k r™+ gjk 6? = Ti k r™+g ik 5J 



cm 



The equations d k gij=digj k = djgi k , d k Y s ^ = diY s - k = djY s ik mean that the components of the metric tensor 
(in these coordinates a k ) are second derivatives of a single function g(&) with respect to the corresponding 
coordinates a k (i.e. gi k = d 2 g/da l da k ), while the components of the affine connection T l - k are also second 

derivatives of a single vector function T(a) (i.e. T s ik = d 2 T s /da t da k , where V s are components of the vector 
function T(a)). Following [5], one can split the affine connection on two parts, i.e. 

V s — n sm n J- f s 

1 ik ~ 2^ ^ mik ^ Jik' 

where the first block \g sm g m ik is a Levi-Civita connection, while the difference of connections /j fe is 
a (l,2)-tensor. Then the equations Y\^g sk =Y s - k g ai =Y\ k g s ^ (see (44)) imply that the tensor f\j k = f^g sk is 
totally symmetric, defining the centroaffine cubic form C = j \jkdx' 'dx^ 'dx k of the hypersurface M N_1 which 
together with the centroaffine metric M = g i jdx l dx : > (satisfying compatibility conditions (44)) uniquely 
characterize a hypersurface. The rest of system (44) are precisely the normalized oriented associativity 
equations, where u 1 = g and all other u k = T k . 

The inverse construction: Now consider an additional variable a 1 and assume that the indices i, j, k 
run from 1 up to N. Since the vector function c(a) has the components determined from (34), we can 
define the quantities ^(a 1 , . . . , a N ) 

c ik = ^fc' i,k = l,...,N, Cj k = Tj k , Cj k = gik, i,j,k = 2,...,N 

and h = (Ar, dr/da 1 , dr/da N ) T . In such a case, linear system (43) is replaced again by more general 
linear system (26). 

Note that these formulas bear considerable resemblance with the formulas for A k for the associativity 
equations (the famous WDVV equation) at p. 36 in [5]. In fact, our formulas reduce to those of Ferapontov 
if the metric g^ is constant in the coordinates a 1 : g^ = fj^ = const. 

Examples: The oriented associativity equations in three-dimensional case with the unity condition 
(33) are nothing but a system of quadratic equations (see (38) and below) 

u b b = v bc w hh - v hb w bc + w 2 bc - w bb w cc , u bc = v cc w bb - v bc w bc , u cc = v\ c - v bb v cc + v cc w bc - v bc w cc . 

Three distinguished versions of the WDVV associativity equations can be singled out by the special choices 
of the metrics g lk (see (6)). If t] u = rj 22 = rj 33 = 1 (and all other entries of if k are equal to zero), then (see 
(5) and (34)) 

„ a 3 b 2 + c 2 „ , 
F = — + ——a + z(b,c), 
6 2 

where the three above quadratic equations reduce to a single one (v = z b ,w = z c and u = (b 2 + c 2 )/2) 

Z bbc Z bcc = 1 + z bbb z bcc ~\~ z ccc z bbc] 



if r] n = 7] 23 = r] 32 = 1, then 



a 3 



F = h abc + z(b, c) 

6 



where the three aforementioned quadratic equations reduce to a single one (u = bc,v = z c , w = z b ) 

1 z ccc z bbb z bcc z bbcj 
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if r] 13 = rf 2 = rf 1 = 1, then 

F = ^{a 2 c + ab 2 ) + z(b,c), 
where the three corresponding quadratic equations reduce to a single one (u — z c ,v — z b , w = b 2 /2) 

Zccc ^bbc ^bbb^bcc- 

The first two examples are associated with the hypersurfaces endowed with flat centroaffine metrics 
(see [5] for details), while the third example is related to a non-flat centroaffine metric, because gbb = 
Ubb, Qbc = Ubc, g C c = u cc and all components of the Riemann curvature tensor do not vanish. 

Thus, we established a link (in fact, an equivalence) between the oriented associativity equations with 
unity and centroaffine geometry in a general non-flat case (cf. [5]). 



8 Inverse Construction 

In the preceding sections we constructed the transformation from symmetry consistent conjugate curvilin- 
ear coordinate nets to the oriented associativity equations. In this section, we briefly discuss the inverse 
transformation. 

• Any solution of oriented associativity equations (1) is associated with the corresponding hydrodynamic- 
type systems (4). 

• Suppose that all characteristic velocities w[ fc )(a) of each hydrodynamic-type system (4) are distinct, 
i.e. the algebraic equations (for each fixed index k) 

det |c} fc -^ (fc) 5:;.| =0 (45) 

have just simple roots. In this paper we restrict our consideration to this semi-simple case only. 

• If an iV-component hydrodynamic-type system has pairwise distinct characteristic velocities (see 
(45)), iV conservation laws (see (17)) and all components of the Haantjes tensor (see [7]) vanish, 
then this hydrodynamic-type system is semi-Hamiltonian (see [17]) and can be written in diagonal 
form (i.e., N Riemann invariants exist). The Nijenhuis tensor for the hydrodynamic-type system 
(23) reads (below in this Section d k = d/da k ) 

N} fc = vp p vl - vld p v) - vlidjvl - d k v% 

and the Haantjes tensor has the form 

For each system (4) with the time variable t s we readily find that the Nijenhuis tensor reduces to 
the form (recall that the summation is over the pairs of oppositely located repeated indices only) 

TVp _ r 1 J _ r 1 J 

1 ^(s)jk ^js^qks ^ks^qjsi 

while the Haantjes tensor becomes 
H* = (P- r m (r q r { —r q r { )-\-r i r m (r q c p —r q c 9 "l-l-r* r m (r q —r q r 9 )-\-r i r p (r q r n —r q r n ) 

(s)jk u js^ks\^ps qms ms^qps J'^ps^ksK^ms^qjs ^js^qms/'^ps^jsK^ksqms ^ms^qks/' ps^nsK^js^qks ks qjsJ ' 
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This expression formally contains eight blocks. However, the fourth and fifth blocks coincide. Thus, 
the Haantjes tensor reduces to the six-block form 



xqi = J> r m ( q i _ q i \ , i q ( m J> _ m J> \ , i r> 



L (s)jk L 'js^ks\ L "ps'-'qms u ms u qps) ~T~ { -'ps L 'ms\ L 'ks { -'qjs ^js^qks) ^ps^nsK^js^qks ^ks^qjs) " 

Then one can rewrite the Haantjes tensor in the following form 

xri r m f) (r 1 r q — r l r q \ 4- r l r m F) (r q — r q r p \ 4- r l r m f) (r p r q — r 9 r q \ 

(s)jk ^js Lj ks u s\^qm'-'ps ^qp^ms) ~ ^ps^ks u s\'- J ms^qj ^mj^qs/ ~ ^ps u js u s\^qs^mk ^qk^ms) 

a_(J P m _ J r m \r p r q 4- (r { r m — r l r m \r p r q 4- r l ( r m — r m \r q 
'K^qm^js ^ms^qj J^ks^pss ' K^ms^qk ^qm^ksJ^js^pss ~ ^qsV^ks^mj ^js^mkJ^pssi 

where each bracket vanishes by virtue of oriented associativity equations (1). So, indeed, the family 
of hydro dynamic-type systems (4) is semi-Hamiltonian. 

In such a case (see (22)), the Riemann invariants can be found by quadratures 

iV 



£ f v\ m) {*)dcr : (46) 

m=l ^ 



where the characteristic velocities w^(a) also (see (45)) satisfy quadratic relations (cf. [4]) 



N 



s=l 



This means that the functions a n (upon inverting the point transformation (46)) can be expressed 
via the Riemann invariants r k . 

Then following Tsarev's construction (see [17] for details), one can compute the so-called Lame 
coefficients (the expressions on the r.h.s. are equal for all values of s) 

d k \nHi = , d fc 4o( a ) / 47 n 
«?.)(»)-«(.) (a) 

Then following Darboux (see [1] for details), one can find the rotation coefficients 

0* = ^. i*k, (48) 

which satisfy 2 system (7), (8) describing symmetry consistent conjugate curvilinear coordinate nets. 
Indeed, since the rotation coefficients are the same for all commuting flows (4), without loss of 
generality consider just N — 1 commuting hydrodynamic-type systems (31). Then (46) reduces to 
the form 

N „ 



1 + E / *{m)(s)da" 

m=2 ^ 



where f( m )(a) depend on the variables a 2 , a N only. Thus, Sv l , m ^(a) = 0, where the shift symmetry 
operator (see (35)) 



d y^dr m d d 

m=l m—1 



2 Since the Lame coefficients Hi are determined up to multiplication by arbitary functions /^(r 1 ) (see (47)), they should 
be fixed by the restriction 5(3 ik = 0. 
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Since the Lame coefficients Hn^ can be found from (cf. (47)) 



dkV % , s \(a) 

we have d k 5\n.H(i)i = 0, i.e. ^lnif^j = Xiij 1 )^ where Xiij" 1 ) are arbitrary functions. Thus, one can 
choose H(i)i = Xi(r l )H(i)i such that Xj( r *) — ®% l n Xi( r *)- Then SH^ = 0, i.e. the corresponding rota- 
tion coefficients (3 ik depend on differences of the Riemann invariants, because 5(3 ik = 0, which follows 
from (48). So, we conclude that any solution of the oriented associativity equations (1) gives rise to 
some solution of system (7), (8) describing symmetry consistent conjugate curvilinear coordinate nets. 



9 The Widest Class of Semi-Hamiltonian Hydrodynamic-Type 
Systems 

This section is devoted to description and integrability of the widest class of semi-Hamiltonian hydrody- 
namic-type systems (27). A general solution of oriented associativity equations (1) leads (see Section 5) 
to a general solution of system (7), (8) describing symmetry consistent conjugate curvilinear coordinate 
nets as well as to the basic set of solutions H^ k (see (11)). In this section we construct general solutions 
of linear spectral problems (9) and (10) which are important for various applications in the theory of semi- 
Hamiltonian hydro dynamic-type systems, for instance, in the generalized hodograph method (see [17]). 
Recall a few major formulas from Sections 2, 3, 4. 

Suppose that the basic set of solutions Hn\ k (see (11)) of the linear spectral problem (see (9) and (10)) 

SH^XHi, d i H h = p ik H i , i^k (49) 

is found for a given set of rotation coefficients (3 ik depending only on differences of the Riemann invariants 
r n (see (8)) and satisfying (7). Then, the basic set of solutions ip^ of the adjoint linear problem (see (9) 
and (10)) is given by (13), where g sn is a non-degenerate symmetric metric which is inverse to (12). 

Our goal is to find a general solution of the above linear spectral problem as well as a general solution 
of the adjoint linear spectral problem 

Hi = di^k = Pki^i: i ^ k, (50) 

and general solutions of both other linear spectral problems (26) and (28). 

Main result: N infinite series of solutions if - s ' , ip'f' (see (11)) can be found in quadratures. 
Consider the expansions in A of solutions for both linear spectral problems (49) and (50), i.e. (see (11)) 

Hi = H?> k) + Afl-pJ + A 2 tff fc) + ^ = 4°> k) + A# fc) + A 2 # fc) + k = 1, N, 

as well as the corresponding expansions of solutions of two other linear spectral problems (26) and (28), 
i.e. 

h = A" 1 ^ + a k + Xc k + X 2 h {2 ' k) + A 3 /i (3 ' fe) + b s = X^Sl + c s k + Xb s (hk) + X 2 b\ m + A 3 ^ 3;fc) + 
where <9j/i (n,fc) = ipf 1 '^ H^yi and dib* n ^ = ip^ H- n ' k \ while (see Section 4) 

N N 

n — °(n,k) ~ / n m 

m=l m=l 



16 



Instead of complicated integration of infinite sets of recursion relations (11) (see (49) and (50)), one can 
easily integrate the other infinite sets of recursion relations (see (26) and (28)), 



N N N N 



dh ^) = J2J2ct m h^ k) da m , db\ n+hk) = Y,Y, c ™ b Uk) dam > i = l,-,N, n = 0,l,... (51) 

m=l s=l m=ls=l 

The construction of all these four recursion relations includes the following steps: 

• Since ipf"'^ = = c^ipf^ (see the proof of tl 
expressions c k can be found in quadratures, i.e 

TV 

dc k s = J2i>m H (s)mdr 



Since ipf"'^ = ip^ = c^ipf^ (see the proof of the theorem in Section 3), where diC k = tjj\ k ^ H^y, the 



m=l 



Then (recall that ^a s = ijj'f'H^i, diC s = ijjf 

N N „ N N / N 

4' k) = E^E / ^ )H (s)mdr m , dc k = £J>« E / ^ )H (s)mdr m 

s=l m=l p=l s=l \m=l 



Thus, we had a basic set of solutions H^ k . Then we reconstructed the adjoint basic set of solutions 
tp^m • Then we found N conservation law densities a k and N conservation law densities c k . So, we 
have found all structure constants c* fc with unity condition (33). 

• Then step by step (see (51)) we can find N infinite series of higher conservation law densities h^ n ' k ^ 
as well as N infinite series of higher conservation law fluxes b^ n k ^ . 

• Since dih^ n ' k ^ = ijj[ n ' k ^H^i and dib s ^ nk ^ = ip[ s 'H^ n ' k \ we can find higher solutions Hf' k \ ip^ ,k \ 

Moreover, we want to find general solutions of linear systems (9) without extra conditions (10). This is 
essential for application of the generalized hodograph method (see [17]). However, just in some very special 
cases general solutions of linear systems (9) can be found explicitly. Nevertheless, using the approach 
suggested in [17], one can see that any initial data for the linear system (9) can be approximated by linear 
combinations infinitely many particular solutions H\ n,k \ tp^ n ' k ^ with appropriately chosen coefficients £ n fc . 

This means that the general solution of an arbitrary semi-Hamiltonian hydrodynamic-type system 
(27) whose rotation coefficients depend on differences of Riemann invariants only can be found using the 
generalized hodograph method (see [17]) and has the form 

oo N 

(n,k) 



xHi + tHi = EE^* F < 



n=0fc=l 



10 Conclusion 



In this paper, we considered three distinguished objects: oriented associativity equations, symmetry con- 
sistent conjugate curvilinear coordinate nets and semi-Hamiltonian hydrodynamic-type systems whose 
rotation coefficients depend on differences of the Riemann invariants only. We have shown that these 
objects are closely related, and thus the knowledge about one of them implies the knowledge about the 
others. 
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